Surface oscillations in channeled snow flows. 
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Abstract 

An experimental device has been built to measure velocity profiles and friction laws in 
channeled snow flows. The measurements show that the velocity depends linearly on the vertical 
position in the flow and that the friction coefficient is a first-order polynomial in velocity (u) 
and thickness (h) of the fiow. In all flows, oscillations on the surface of the fiow were observed 
throughout the channel and measured at the location of the probes. The experimental results 
are confronted with a shallow water approach. Using a Saint- Venant modeling, we show that the 
flow is effectively uniform in the streamwise direction at the measurement location. We show 
that the surface oscillations produced by the Archimedes 's screw at the top of the channel persist 
throughout the whole length of the channel and are the source of the measured oscillations. This 
last result provides good validation of the description of such channeled snow flows by a Saint- 
Venant modeling. 
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1 Introduction 



In the past, many models have been used to describe the behavior of dense snow avalanches 
fVoellmv, 1955:'Perla ct al.l. ll980l : lDent and Lanel . ll983t[Norem et aLl . ll987t[Savage and Huttei . 
[l989;.Naaim et al. . 2003). Since snow avalanches are very destructive and cannot be controlled, 
it is very difficult to fully t est these models experimentally. Direct measurements on avalanches 
do exist ( Dent et al.l . ll998l ) but are rare and must be comple ted by " laboratory expe riments" to 



study the rheology of snowl Nishimura and Maeno , 1989, : .Casassa et al. , 1991 ; Kern et al. . 200^ . 



The main problem is that snow grains are very fragile. This brittleness makes it very difficult to 
use a lar ge amount of sn ow in the l aboratory. In consequence, as described in previous papers 
(jBouchet . .2003a; Bouchet et al.l . l2003h) we chose to work with an in situ rheometer, to use snow 



directly from the field. To stay as close as possible to naturally occurring flows in avalanches, 
we studied free surface snow flows in an inclined channel. 
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Experiments over a wide range of slopes and thicknesses were perf ormed to measure vel ocity 



profiles and friction laws in the flow. Results arc fully reported in fBou chet et al.1 . 120041 ) and 



are briefly presented in section [2721 One of the main purpose of the present article is to test the 
obtained snow flows with regards to a commonly used Saint- Venant model to see if the model 
can help in understanding and predicting the behavior of such flows. 

As described later in the article the snow flow consists of a short flrst phase (2 to 3 seconds) of 
transient regime followed by a fully developed regime which lasts until the end of the experiment 
(around 30 seconds). The fully developed part of the flow is made on the one hand of a 
steady component non uniform in the streamwise direction but which tends to become uniform 
asymptotically (x-uniform asymptotic regime) as will be shown in the article. On the other 
hand, added to the steady component is a periodic one which will be studied in detail when 
looking to the surface oscillations that are present on top of the flow. The whole experimental 
and numerical results presented in the present paper are all focused on the fully developed 
regime leaving apart the transient one. The steady component of the fully developed flow 
is studied in detail in this article using the analytical and numerical results. As mentioned 
before, in all cases, oscillations on the surface of the flow were observed. The question of the 
origin of these oscillations has not been addressed up to now and is another main concern of 
this paper. Are t hey ro ll waves intrinsic to the flow like those observed in granular flows by 
Louge and Keasj ( 2001 )? Are they engendered by the feeding system as direct observations 



suggest? Numerical simulations and analytic solutions are used to get a better understanding 
of the flow characteristics. 

Section [2] briefly presents the experimental set-up and results. The simulation technique is 
presented in section [31 it is based on the resolution of the shallow water Saint- Venant equations. 
The behavior of the flow with increasing streamwise distance, and in particular the way it 
approaches the x-uniform, fully developed regime, is then studied numerically and analytically. 
An estimate of the distance needed to attain the asymptotic regime is derived. Section|4]concerns 
the study of the oscillations on the surface of the flow, and the validation of the model. 



2 Experimental set-up and results 

The experimental facility and results are summarized in iBouchet et al. (I2nn3bl . l2004l ). A more 
detailed technical description of the set-up, and complete results, can be found in Bouchel] 
( 2003a[) . 



2.1 The experimental device 

A photograph of the experimental apparatus is shown in figure [TJ The device consists of a 
10 m long, 20 cm wide, 20 cm high channel, fed with snow by an Archimedes's screw from a 
hopper at the top of the channel. The screw is driven by a thermo-hydraulic engine, so that its 
rotation speed / is constant during an experimental run. Such a system is able to produce fully 
developed flows with a typical duration of 30 seconds. The floor of the channel is covered with 
sandpaper to create a rough surface. The slope angle of the channel, 9, ranges from 27° to 45°, 
in steps of 2°. A second adjustable parameter is the mean flow thickness {h) which is controlled 
by the value of the rotation speed f (f ranges from rpm to 60 rpm) . Twelve different flows were 
studied. Many types of snow can be found in the field. As a first study the repeatability of the 
experiments was privileged and thus flows were performed in the following conditions. In each 
case: 

- to avoid snow melting the temperature of both snow and air was below — 15°C; 

- snow grains were small and rounded with diameters between 0.2 and 0.4 mm; 

- snow was sieved with a 3 cm mesh. 
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Figure 1: The experimental set-up is designed to perform free surface snow flows in a channel 
(length: 10 m, width: 20 cm, thickness: 20 cm). The coordinate system associated to describe the 
flows (see figure H]) is the following: the x axis is along the main direction of the channel, orientated 
downward; the z axis is perpendicular to the bottom of the channel, orientated towards the surface 
of the flow. 
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The velocity profile was measured 6 m downstream from the top, on one of the sidewalk. The 
flow thickness was measure d at the same loca tion {Hl) and also 1.5m upstream (/i^-i s). 

Velocities are deduced ( Dent et al. . 1998[ ) from the delay between signals at two identical 
opto-electronic sensors, one 2 cm downstream from the other. Each of these sensors consists of 
a photo-transistor and an LED. The latter emits an IR light beam which is reflected by the 
snow and collected by the photo-transistor. The signal produced by the photo-transistor is thus 
characteristic of the snow pack flowing in front of the sensor. Because the distance d between 
the two sensors is small enough to avoid deformation of the snow pack, the signals are identical 
but delayed by an amount of time r. A cross-correlation method enables measurement of r, and 
hence to deduce the velocity v ~ d/r. 

The flow thickness measuring device (Leuze ODSM/V-5010-600-421) consists of an optical 
distance sensor positioned above the flow. It emits an IR light beam (a few degrees) which is 
reflected by the free surface of the flow and collected by a CCD sensor. The measured position 
of the reflected beam on the CCD sensor, is a function of the distance between the sensor and 
the surface. Knowing the distance between the sensor and the bottom of the flow allows the 
calculation of the flow thickness. 

The main characteristics (slope angle, average thickness and frequency of the thickness os- 
cillations) of the flows are listed in table [TJ The frequency F of the oscillations was obtained by 
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F: frequency of the 
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thickness oscillations (Hz) 
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37 
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0.42 


39 
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0.38 


39 
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0.50 


39 


9.35 


0.68 



Table 1: Slope angle, average thickness and frequency of the thickness oscillations for each analyzed 
flow. 



a Fourier transform of the signal received from the thickness's sensor. 

2.2 Thickness oscillations, velocity profile and mean velocity 

Initially the channel is empty. The beginning of the experiment consists of a transient phase, 
while the front of the flow goes throughout the channel. After this transient phase, as shown in 
figure [21 at any given location the flow thickness oscillates in time around a mean value < hx > 
defined as a time average (steady component of the flow). As also illustrated in figure [21 the 
variations of follows those of /i;-i.5 with a constant delay r: h^lt) « /ij;,_i.5(t — r). It thus 
appears that the steady component of the flow is x-uniform at the location of the probes. The 
frequency of the thickness oscillations mentioned above is always close to the rotation speed of 
the screw, as illustrated by figure [3] Indeed, the rotation speed of the screw can be calculated 
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Time (s) 



Figure 2: Example of thickness variations. Hl (full line) is measured 1,5 m downstream from hz-i.b 
(dashed line) . The figure corresponds to a laps of 10 seconds taken from a fiow lasting approximately 
30 seconds. 




Figure 3: Frequency (F) of the oscillations observed as a function of the screw rotation speed (/). 
The straight line is given hy F = f. 
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according to the formula: / = V^^ x Q with Q the flow rate (obtained from the velocity profile 
as described below) and V « 0.14 m"^ (information given by the manufacturer). / is calculated 
with a 10% error. 

Special attention was payed to side wall effects on the flow. First the existence of a transverse 
velocity gradient was qualitatively looked for. To that purpose the flow was filmed from above. 
The observation show a di splacement of the snow as a block without slowdown toward the 
walls. More guantitativelv iLouge and KeastI ( 200l[ ) measured velocity gradients in a granular 



channeled flow and obtained less than 10% in velocity between the center of the flow and the side 
walls. Because in the present case the material is cohesive the gradient should be eve n less than 
these 10%. Side walls could also be responsible for a global slowdown of the flow. Erickssoiil 



performed measurments of friction coefficients between snow and different materials. 



The friction coefficient between snow and plexiglas was measured to be around 0.05 which is 
ten times lower than the friction between snow and the sandpaper sticked on the bottom of the 
channel. Thus, it can be assumed that the effect on the mean velocity is negligible. From these 
two points, the ffow can be considered as uniform in the direction perpendicular to the flow. 

The mean velocity profile is written V{z), where z is the distance from the bottom of the 
channel. The r esults for the twelve different flows are all satisfyingly fitted by a linear profile 
(|Bouchet et all . [2004) : 



V{z-e,h) = V,+Toz (1) 

where: 

Vs (9) = Vo [tan 6 - tan 6*0] (2) 

and 

Vo =9.1 m.s-i 

tan 6*0 = 0.28(6*0 = 15.6°) (3) 
Fo = 12.5 s-i 



Such a linear profile i s not us u al but has already been o bserved in some granular flows (jPouadv et al 



2002t iBonamv et al.L [2002a[ ). iBonamv etHl (|2002bh explained this uniformity of the velocity 



gradient by the presence, in the flow, of aggregates of many different sizes wi thout one typi- 
cal scale. In the present flow, aggregates were also detected (|Bouchetl . l2003al) . It should be 



mentioned that the behavior of the flow near the ground was not explored. Indeed, the flrst 
transducer is located more than ten to twenty particles diameters above the ground. Thus, it is 
here impossible to know if there is a basal sliding or a strong basal shear. The presence of the 
sandpaper on the bottom makes the second hypothesis the more believable. Nevertheless in the 
present article we are only interested in the averaged velocity over the flow thickness which is 
quite insensitive to this basal behavior. Moreover, during the numerical calculation the shape 
factor is taken equal to one (see section [3]). 

The integral of equation ([1]) gives the averaged velocity over the flow thickness: 



or, using equation [2l 



u{9,h)^V{e) + ^roh (4) 



u = Vo [tan e - tan 6*0] + ^^oh. (5) 



2.3 Friction law 

Because the mean flow near the probes is fully developed and x-uniform, the friction force exerted 
by the channel on the flow balances the weight of the snow. The effective friction coefficient /i, 
defined as the ratio between the shear stress and the normal stress at the bottom of the channel, 
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Figure 4: Sketch of the flow and of the appHed forces. 



is thus equal to tlie tangent of the slope angle 0. Using equation lU, the friction law can be 
written as: 

u h , . 

M = Mo + - (6) 

with 

^0 = tan 6*0 0.28 , , 

Ho = 2Fo/ro = 1.45 m ^ ' 

It must borne in mind that this friction law is here established for the particular dynamical 
conditions of a fully developed and x-uniform flow. Another point is that this friction law has 
been established in a given range of slope angles and flow thicknesses corresponding to those 
of the experiments and should not be extrapolated to some others too far as for example too 
large values of the thickness. The non-transient phase of the following simulations is performed 
close to this range of validity. The initial transient phase corresponds to small flow thicknesses. 
Thus, the friction law is supposed valid even during this phase in our simulations. 

3 Study of the steady component of the flow 

In this section we study in detail the behavior of the steady part of the fully developed flow. To 
that purpose we use the laterally unconfined Saint- Venant equations solving them first numeri- 
cally and in a second time analytically. 

3.1 The Saint- Venant equations 

Originally used for hydraulic stu dies of shallow water flow s, the Saint- Venant equations were 
first applied to granular flows bv ISavage and Hutter (Il989l ). This approach is now commonly 



employed to model dense snow flows. The equations are derived by averaging the mass and 
momentum conservation equations over the flow thickness. Two assumptions are employed: 

• the typical flow thickness is assumed to be small compared to the streamwise variations 
(short wavelength hypothesis); 

• the flow density is taken as constant and uniform. 



7 



Hence, the dominant component of flow velocity, u, is parallel to the bottom of the channel. 
The flow is described by its thickness h{x,t) and average velocity u{x,t) = {^^ u{x, z,t)dz) /h. 
Mass and momentum equations for a thin slice can be written: 



dh 

at 



dhu 



dhu 
dx 

dhv? 
I 

dx 



gh cos 9{tan9 — /i 



"dx' 



(8) 
(9) 



where 



• a is the shape factor: a = h u^dzj 

• fc is the ratio between the normal stress parallel to the ground, a^x , and the normal stress 
perpendicular to the ground, azz- 

In what follows, the stress tensor is supposed to be isotropic and thus k is taken equal to 
one. In granular flows, the precise value of k does not appreciably affect the flow behavior. For 
commonly encountered velocity profiles (uniform, linear, parabolic), a takes values between 1 
and 2. The flow is also insensitive to the precise value used for a; here we take a = 1. Finally 
friction law ([6]) is used. 



3.2 Description of the numerical solution procedure and result 

The 6-meter long channel upstream of the measuring device is discretized in a; with points every 
five centimeters. The time step ia dt — 5 x 10~^ s, and we run the model for thirty seconds. 
The numerical scheme used is an explicit finite difference MacCormack scheme. 

Both steady-inlet and initial conditions are needed for h and u. As regards the steady-inlet 
conditions {h{0), u{0)), once the fiow becomes steady equation ([8]) imphes a constant volume 
fiux: 

h{0)u{0) = h{L, t)u{L, t) = Q (10) 

Experimental values of and ul (given by ([5])) allow Q to be determined. h{0) is taken as an 
input parameter for the simulation and m(0) calculated using m(0) = Q/h{0). These constant 
values of h{0) and u{0), which apply strictly to the steady component of the flow, are imposed 
throughout the calculation. As can be observed in section p.3|) the value used for h{Q) has no 
influence on the asymptotic flow of the end of the channel. 

The values used for 9 and are the mean of those encountered during the experiments: 

/iL = 7.6 cm (11) 

and 

= 35° (12) 

When the model is run, an initial unsteady, non-uniform phase is observed, corresponding to 
the passage of the front along the channel (see flgure[5]). This phase is not discussed because the 
Saint- Venant equations cannot deal well with rapid changes in thickness and velocity. The only 
thing that can be noticed is that the passage of the transient phase lasts around 2 to 3 seconds 
both in the experiments and in the numerics. Once the front has passed, a fully developed 
steady regime is observed (see for example figure [5]) . Near the end of the channel the flow is 
found to be x-uniform. This is an important point for the exploitation of the experimental 
results in general and of the friction law in particular. This regime is examined in more detail 
in the following section. Just for information, the initial values are h(x > 0, t = 0) = 1 cm and 
u{x > 0, t) = 0. The non-zero value of h was introduced to avoid a division by zero. We choose 
these values for the sake of simplicity and because the fully developed flow (the one which is 
relevant to study here) is not sensitive to them. 
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Figure 5: Saint- Venant model's results using ho = 15cm. thickness h (in m), velocity u/W (in 
m/s) and flow rate uh (in m^/s) are plotted versus position x. These results concern the beginning 
of the flow when the front is crossing the channel before the establishment of the fully developed 
regime, (a): t = 0.5 s, (b): t = 1 s, (c): t = 1.5 s. 
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Figure 6: Saint- Venant model's results using ho = 15 cm, thickness h (in m), velocity n/10 (in m/s) 
and flow rate uh (in m^/s) are plotted versus position x{m). These results are obtained in the fully 
developed regime. 



3.3 Analytical study of the steady component of the flow 

In this section, we solve analytically the Saint- Venant equation for the steady component of the 
flow. This enables us to quantify when the asymptotic x-uniform regime is reached. For steady 
flows, the Saint- Venant equations (equations ^ and (O) take the form: 



uh ^ Q (13) 

dhu^ u h dh 

a — - — = (7acos6'(tant^ — /Lto — — + — k—) (14) 

dx Vq Hq dx 

3.3.1 Values of the x-uniform asymptotic regime 

Let hoc and Uoo be the flow thickness and flow velocity when the x-uniform, asymptotic regime 
is reached. From equations (jl3p and (|14p . h^o and Mqo are given by: 

hocUo^^ Q 

Uoo hoc ^ a (15) 

/"0 + 77 — = tan6' 

VO -no 

This system admits only one set of positive solutions: 
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_ffo(tan0 — fio) 



yo(tan6' - ^o) 



(-1 + 



(16) 



;i + VTTi) 



where 



e = 



4Q 



VoHo{t&n0 - fio)^ 

3.3.2 Distance needed to reach the x-uniform asymptotic regime 



(17) 



We next study the development of the flow with streamwise distance x. To that purpose 
we integrate the differential equations (fT5|) and (fH|) analytically using the following reduced 
variables: h = (h— hoo)/hao, u — (w — Uoo)/uoo, with x — (x + xq) / hoa with xq a constant value. 
Details can be found in appendix 1. The integration leads to 



with 









k — aFr'^ 




In 




h 




h + i 



Q_Fr^(l-l//3^) 







h+l + f3 





-aFr 



2I + /3 1 



P /..I 



(18) 



and 



Xq 



Pj.2 

f3 



l + (3 



x{ln[\h{Q)\ 



k — aFr' 



4Q 



VqHq (tan 61 - ^o) 



ghoc cos 9 

UqqHq ^ (1 + yi + e)^ 



(19) 



>(0)+ir^'^'(i-i/^')|/i(0)+l+/3|^'5+aFrV^=]„^£l!(i±i^^ 
' ^ ^ ' I w HI J ^(;j(o) + l) ' 



(20) 



Using equation (|18p we plot /i as a function of x in figure[7l The entire curve has three branches. 
Because the thickness of the flow is limited by the vertical size of the channel {h < 20 cm) 
and because the flow thickness decreases experimentally throughout the flow we only keep the 
part of the curve which fits these conditions. Looking at the curve, one can see that when x 
becomes large enough, the flow reaches asymptotically its x-uniform state. This confirms what 
the numerics seem to predict. The question is now to know whether the asymptotic x-uniform 
regime is reached before the end of the channel. To that purpose we simplify equation p8|) for 
the case of large x. It can be rewritten in: 



(1 + /3) 



'exp 



(l+/3)fcoo 

Ho{aFr^ - k)' 



(21) 
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Figure 7: Plot of the relative deviation of the flow thickness /i as a function of the relative distance 

X. 



The typical length scale is then: 



L = -, = 4 m 

(1 + /3) 



(22) 



when taking the experimental value Q = 0.32 m^.s^^. It is therefore relevant to consider the 
flow as x-uniform at the location of the measuring device, 6 m downstream from the top of the 
channel. 



4 Origin of the oscillations - validation of the model 

In this section we show that the model gives an explanation about the origin of the oscillations 
observed in the experiments. 

4.1 Frequency of the oscillations 

The oscillations induced by the screw which feeds snow into the channel are modeled by the 
boundary condition h{0,t) = ho + /ii cos(wt) with / = cj/27r the frequency of the screw. To 
study what happens at the end of the channel we use the following scheme: 

• First, the model is run for 5 seconds with h{0,t) = Hq + hicoscut to allow time for the 
initial transient phase to have finished. 

• After these 5 seconds the model is run for a further 30 seconds and the temporal Fourier 
transform of h{L,t) is taken. 

An example of the results is shown in figure[51 which also gives the result for h{0, t). In addition 
to the fundamental frequency f and the zero frequency peak representing the mean value, the 
harmonics 2/, 3/, ... are also present in h{L,t) due to nonlinearity of equations ([5]) and 
The amplitude of these frequency components is always less than 10% of the amplitude of the 
fundamental peak. 
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Figure 8: Fourier-transform amplitude of the thickness of the flow obtained by the Saint- Venant 
model. Fourier-transform amplitude of h{0, t), Fourier-transform amplitude of h{L, t). 
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Figure 9: Fourier-transform amplitude of h{L,t) with h{L,t) measured during one experiment. 



4.2 Comparison between model and experiments 

The Fourier transform of h(L,t) in the experiments also shows the presence of the frequencies 
2/, 3/... (cf figure ([9])). However, the ratio of the amplitude of these secondary peaks to that 
of the fundamental differs in the model and in the experiments. In both cases (numerical and 
experimental) the secondary peaks are weak and can thus be neglected. 

Let us now compare the ratio of the amplitude of the fundamental peak to that at zero 
frequency. This ratio is plotted in figure \TU[ With ho = 15 cm and hi = 4 cm the model gives 
good results compared with the experimental values. The curves given by the extreme values 
of the parameters are also plotted on figure [TO] {ho — 12 cm, hi = 6 cm (eg) h^/hi = 2 and 
ft-o = 16 cm, hi — 2 cm (eg) h^/hi — 8). By "extreme values" we mean the larger or smaller 
values of the parameters the experimental observations allow to use. When plotting the curves 
what appears to be important is not the values of Hq and hi separately but the ratio ho/ hi. 
Indeed, with the same ratio, two sets of values (/loa, ^la), [hot, hn,) give merely the same 
curves. This is the reason why we have plotted only two "extreme" cases. One can see that if 
the experimental values are close to the predictive curve they are also all included in the space 
between the two extreme curves. This allows for different conclusions. The majority of points 
(those that are close to the predictive curve) show that the model well predicts the behavior of 
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Figure 10: Ratio between the amplitude of the fundamental peak and the amplitude of the peak 

at zero frequency versus frequency. * results from experiments, results of the Saint- Venant 

model with ho = 15 cm and hi = A cm, results of the Saint- Venant model with /iq = 12 cm 

and hi = Q cm, results of the Saint- Venant model with Hq = 16 cm and hi = 2 cm. 

the flow. The model is run in each case with fixed parameters {0 = 35°, Hl = 7.6 cm). These 
values are the mean of the 11 first experimental values (see table [T]), the twelfth amplitude 
was not usable. Hence it seems reasonable that some of the points are not strictly close to the 
predictive curve. Nevertheless they are still between the extreme curves. All of this confirm the 
hypothesis that the model well predicts the phenomenon and the flow and that the oscillations 
are due to the screw. 

5 Conclusion 

The dense snow flows studied during the experiments exhibit thickness's oscillations. Though 
experimental facts suggest that they are generated by the feeding system of the set-up, the origin 
of these oscillations remain a priori uncertain. The friction law for dense snow flows deduced 
from the experiments is used in numerical resolutions of Saint- Venant 's equations. The initial 
conditions given to the thickness at the entrance of the channel consist in sinusoidal oscillations 
with an amplitude hi, around a mean value ho- The simulations performed lead to the following 
results: 

• the mean flow is x-uniform at the location of the measuring device. This is a crucial point 
for the validation of the results of the experiments in which the hypothesis of x-uniformity 
was the central hypothesis; 

• the oscillations engendered at the channel's entrance propagate and are still present with 
the same frequency (harmonics can be neglected) at the measurements' position; 

• numerical and experimental results are in a good agreement when using an oscillation ratio 
hi/ho = 4/15. From direct observations, this ratio seems very reasonable; 

• the relative intensity of these oscillations (given by the FFT of h{x, t)) reduces from about 
15 % at the top of the channel to approximatively 5 % to 10 % at the studied position. 
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The results presented in this paper agree with the experimental observations. They suggest 

that: 

1. in the experiments, the mean value of the thickness of the flow at the entrance of the 
channel can be changed as wanted without affecting the x-uniformity of the flow at the 
end of the channel and thus the vahdity of the measurements; 

2. reasonable oscillations generated by the feeding system can propagate all along the channel 
and have the same intensity as the one observed in the experimental results; 

The present study reveals that, as many other flows, such channeled flows of dry snow can be 
described using a Saint- Venant like model. Such numerical studies can help in the future in 
preparing other configurations of the device to test other caracteristics of those channeled snow 
flows and thus rheology of snow. 
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Appendix 1 

Using the mass conservation equation (equation (|13p ). equation p4p can be rewritten in: 



dx 



kg cos 9 , 
aQu H h- 



gh cos 9 



tan 9 — /iQ 



Vo Ho 



(23) 



Let Su and Sh be the deviations from the x-uniform regime: 



u = Uc 

h — 



Su 
Sh 



Equation becomes: 



d_ 

dx 



aQSu + {2ha,Sh + Sh') 



gihoa + Sh) COS I 



Sh Su 



(24) 



(25) 



We now turn all lengths and velocities without dimension by dividing them respectively by 
/loo and by Moo and define: 



X — xjhoQ 

h = [h - hoo)lhoo 

U = (u-Moo)/Moo 



(26) 



Mass conservation (equation (fT3|) ) implies: 



hu + h + u = Q 



(27) 



At this stage let's notice that if h tends to zero, because of equation (|77|) . u also tends to zero 
(consequence of mass conservation) (eg) the x-uniform regime is reached. Hence in what follows 
we focus our study on h. For its part, equation (|25p becomes: 



d_ 

dx 



aFr'u + kh+ ]^kh' 



= [l + h)[^h-^u) 

no Vo 



(28) 



with Fr the Froude number of the x-uniform regime: 

Fr^ ^ , (29) 

V.g/ioo cos 9 

Substituting, in equation (j28[) . u by its expression obtained from equation (j27p and intro- 
ducing the parameter j3 = UooHol {haoVo), gives: 



d_ 

dx 



kh+-h' - aFr'rJ^ 
2 h+1 



= -^h(h + l + (3 
Ho 



which can be rewritten in: 
h + l 



aFr'-^ 



1 



h{h+l + (3) h{h + I + I3){h + If 

A straightforward integration of equation PT|) leads to: 



dh 



-dx 



(30) 



(31) 
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h+1 
2I + /3 



QFr^(l-l//3^) 



aFr 



h+1 -no 



= ^o(tan'^-Mo) [i + yiTi] 

Foi?o (tan^ - /xq)^ 
ghoo cos 

Mooffp ^ (1 + yi + e)^ 
hooVo e 
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